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Inelastic neutron scattering data from a twinned single-crystal of YBa 2 Cu 3 06 95 are presented 
that show a distinct a-b plane anisotropy in the oxygen vibrations. The Cu-0 bond- stretching type 
phonons are simultaneously observed along the a and b directions due to a 4 meV splitting arising 
from the orthorhombicity. The present results show the bond- stretching branch along b (parallel 
to the chain) has a continuous dispersion, while the branch along a is discontinuous, suggesting a 
possibility of short-range cell-doubling along a. Furthermore, the LO mode along a is split in en- 
ergy from its TO partner at non-zero ^-vectors, while the b mode is not. These results imply 
strong anisotropy in the electronic screening and a one-dimensional character in underlying charge 
fluctuations. 

PACS numbers: 74.25.Kc, 63.20.Kr, 71.30.+h, 74.20.Mn 



In recent years, evidence has been building that the 
lattice may play a non-trivial role in the transport and 
superconducting properties of high-temperature super- 
conductors. In particular, certain modes (the high-energy 
Cu-0 bond-stretching modes) show strong softening 
with doping, and their frequencies are strongly and 
abruptly reduced when crossing from the insulating to 
the metallic phase as a function of doping [1-4]. Photo- 
emission also suggests strong electron-phonon interac- 
tion for these modes [5,6]. The compendium of meas- 
urements seem to support an anisotropic, locally inho- 
mogeneous charge model, which may be related to the 
stripe model [7,8]. In the present Letter, we show that the 
dispersion of Cu-0 bond- stretching phonons in ortho- 
rhombic YBa 2 Cu 3 695 (YBCO) are very different along 
the a and b directions, with the dispersion along a being 
discontinuous halfway to the zone boundary. This is also 
the modulation wave vector direction for incommensurate 
spin fluctuations observed in untwinned single-crystals 
[9] establishing a connection between the Cu-0 bond- 
stretching modes and the strioe structure. 



Inelastic neutron scattering measurements were per- 
formed on a twinned single-crystal of YBa 2 Cu 3 695 
weighing - 100 grams. Measurements were focused on 
the Cu-0 bond-stretching type phonons along the a/b- 
direction in the twinned crystal. Untwinned single- 
crystals of sufficient size for phonon measurements are 
not available. The a-direction is called £ and the b- 
direction A according to standard notation [10]. YBCO 
has a distinct orthorhombicity due to the presence of 
CuO chains running along the ^-direction. In addition, 
the YBCO crystal structure contains two Cu0 2 planes 
separated by 3.38 A. The high-energy bond- stretching 
phonons are oxygen modes polarized in the Cu0 2 plane 
and can have either in-phase displacements between the 
two planes (T\ symmetry) or out-of-phase displacements 
(r 3 symmetry). Previous measurements have shown that 
these phonon modes have strongly reduced frequencies 
compared to the insulating parent compound YBa 2 Cu 3 6 
between q = (x/2a, 0, 0) and (n/a, 0, 0) and the corre- 
sponding wave vectors along b [4] . 
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FIG. 1 . Neutron scattering intensity map by pulsed neutron ine- 
lastic scattering performed at the ISIS using the MAPS spec- 
trometer at T = 110 K. Data are shown as a function of energy 
transfer and momentum transfer (g), in the range 2.5 < q x < 3.5, 
-0.05 < q y < 0.05, in the unit of reciprocal vector, 2n/a. The data 
were divided by q 2 , after subtracting the background, since the 
phonon intensity is proportional to q 2 . The data were taken at T 
= 1 10 K. White lines are due to gaps in the detector coverage. 



Measurements were performed by neutron time-of-flight 
spectroscopy on the MAPS spectrometer at the ISIS facility 
of Rutherford Appleton Laboratory using an incident en- 
ergy of 118 meV. The sample orientation had the a/b-axis 
perpendicular to the scattering plane and the oaxis was 
rotated 41° away from the incident beam direction. MAPS 
data were analyzed using the Mslice program. Additional 
measurements on the same crystal were performed by tri- 
ple-axis spectrocopy on the HB-2 and HB-3 instruments at 
the High-Flux Isotope Reactor at Oak Ridge National Labo- 
ratory. For the triple-axis measurements the final energy 
was fixed at 14.87 meV and the sample-to-detector horizon- 
tal collimation elements were 48'-60'-40'-120'. 

Due to the presence of the Cu-0 chains, some degree of 
anisotropy is expected in the lattice dynamics in the Cu0 2 
plane from orthorhombicity. At the Brillouin zone center, 
Raman scattering on untwinned single-crystals has identi- 
fied that the B 2g (a-axis polarization) and B 3g (Z?-axis) modes 
which connect to the ( r 3 ) bond-stretching modes are split by 
4 meV [11]. We note that the in-plane displacements of the 
chain oxygen modes are not Raman active, except perhaps 
weakly due to disorder, therefore the above mode assign- 
ments must be correct. This splitting is also confirmed by 
neutron scattering on twinned crystals [12] and the results 
presented here. Neutron scattering shows the same splitting 
for the B 2u and B 3u modes, which connect to the T\ modes 
[13]. The T 3 modes around the reciprocal lattice point (3, 0, 
0) (in the unit of reciprocal lattice vectors) are shown in 
Fig. 1 where the B 2g and B 3g modes are observed at 72 and 
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FIG. 2. The dispersion of the Cu-0 bond- stretching phonon 
modes in YBa 2 Cu 3 6 95 at T = 7 K, determined by curve fitting to 
the MAPS data. The shift between the modes seen in Fig. 1 is 
corrected here. The left-hand side describes the LO phonons and 
the right-hand side the TO phonons. The branch that is denoted 
B 3g at the zone-center (o) has polarization along the b axis (paral- 
lel to the chain). This mode is smoothly dispersing and TO (£) 
and LO (A) are nearly degenerate. On the other hand, the B 2g L 
branch that has polarization along a ( ) is discontinuous at the 
middle of the zone, and joins the B 3g -A mode at the zone-edge. 
The LO frequency of the B 2g J L branch is significantly lower than 
the TO (A) frequency for non-zero q. Lines are guide to the eye. 



68 meV, respectively [14]. Note that the dispersion maxi- 
mum for the B 2g mode occurs around (3.02, 0, 0) with the 
provisional lattice constant of 3.85 A, while that for the B 3g 
mode is located at (2.96, 0, 0). This is because a is smaller 
than b by 2 %, confirming that each mode propagates along 
a different orthorhombic axis. The phonon splitting due to 
orthorhombicity is consistent with the higher frequency 
mode having a shorter unit cell distance. 

The splitting of the dispersions of the A and £ modes 
allow to follow each mode independently, in spite of the 
fact that the crystal is twinned. As shown in Fig. 2 the lower 
energy A branch disperses strongly from 68 meV down to 
55 meV near the zone boundary Y-point (0, 7i/b, 0) where it 
mixes with an oxygen bond-bending branch of the same 
symmetry. The E branch, on the other hand, disperses from 
72 meV to 68 meV at the (7d2a, 0, 0) point where the inten- 
sity rapidly disappears. The intensity of the L branch reap- 
pears near the zone boundary at ~ 55 meV. This discontinu- 
ous nature of the bond- stretching phonon branch can arise 
from strong electron-lattice coupling and may be under- 
stood as a dynamic cell-doubling [1,15,16]. Here we show 
that this gapped behavior of the bond- stretching branch 
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occurs only for bond- stretching phonons propagating along 
the a-direction, while the /^-direction has a continuous dis- 
persion. This suggests that the formation of short-ranged 
charge fluctuations with length scale 2a occur only along 
the a-direction. 

The one-dimensionality of the electron-lattice interaction 
is further enforced by examining the transverse modes (with 
T 2 symmetry) emanating from the B 2g and B 3g phonons. The 
behavior of the transverse modes is shown in the right-hand 
side of Fig. 2. The transverse mode along A (with polariza- 
tion along the a-direction) remains dispersionless at 72 
meV until at least (0, 7l/2b, 0), the extent of our measure- 
ment range. The transverse mode along £ has a dispersion 
similar to the longitudinal bond-stretching A mode. As 
expected for a metal, the LO/TO splitting is zero at the zone 
center due to electronic screening. The LO/TO splitting 
remains very small for oxygen displacements along the b- 
direction, however, it increases with q along a. This implies 
that the phonons polarized in the b -direction remain well 
screened and more metallic than those polarized in the a- 
direction. In other words, the electronic dielectric function, 
€ei(q, (d), is different along a and b and the low-frequency 
charge fluctuations in YBCO are very anisotropic. More- 
over, the LO frequency is lower than the TO frequency for 
non-zero-^, which is anomalous in the spirit of the 
Lyddane-Sachs-Teller relationship [17]. The effect of elec- 
tronic screening due to charge fluctuations was investigated 
via the calculation of £&(q, cd) by Tachiki and Takahashi 
[18]. They have shown that overscreening due to the vi- 
bronic state, in which charge and lattice vibrate at the same 
frequency with some phase shift, results in the similar sof- 
tening of the LO mode for q ^ 0. 

The dispersions of the LO bond- stretching phonon 
branch were found to show no apparent temperature de- 
pendence, but the phonon intensity was found to show sig- 
nificant change with temperature, beyond what is expected 
for thermal excitation. Figure 3 shows the difference in 
intensity between 7 K and 110 K. There is a transfer of the 
phonon intensity from a position halfway to the zone 
boundary (in the energy range from 55 - 70 meV) to the 
zone-center (70 - 75 meV) and zone-boundary (50 - 55 
meV). A very similar difference pattern was observed also 
in the triple-axis data [19]. This transfer in spectral weight 
can be expressed in terms of the intensity at q = (fl/2a, 0, 0) 
as the difference in the average intensity between the Range 
1 (51 - 55 meV) and the Range 2 (56 - 68 meV). Fig. 4 
shows the difference, 1(1) - 1(2), obtained with the triple- 
axis spectrometer. It grows below T c , just as the supercon- 
ducting order parameter. 

The anisotropy of the phonon dispersion and implied 
anisotropy of electronic screening suggests that charges 
have low-frequency fluctuations which are strongly anisot- 
ropic. They could possibly originate from an array of ori- 
ented one-dimensional metallic channels separated by less 
metallic regions, as in the stripe model [8,20,21]. Then, 
ionic displacements along the metallic channels are well 
screened, whereas displacements perpendicular to the chan- 
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FIG. 3. Difference in the scattering intensity between 7 K and 
110 K, along h, determined by pulsed neutron scattering with the 
MAPS at the ISIS. To increase the statistics the data for q = 2.5 - 
3.0 and q = 3.0 - 3.5 were averaged after mirror operation. Be- 
low T c the intensity increases in the energy ranges 50 - 55 meV 
and 67 - 75 meV, while decreases in the energy range 55 - 67 
meV. 
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FIG. 4. Temperature dependence of the difference in the neutron 
scattering intensity at q = (3.25, 0, 0) between the energy range 1 
(51-55 meV) and range 2 (56 - 68 meV). The plot shows 1(1) - 
1(2), determined by the triple-axis spectrometer. The supercon- 
ducting transition temperature of this sample is 93 K; the data 
show a striking resemblance to the temperature dependence of 
the superconducting order parameter. 



nels are poorly screened due to an energy gap for holes 
jumping out of the channels. The incommensurate spin 
fluctuations signify the presence of such charge and spin 
modulation [9]. Inelastic neutron scattering measurements 
with untwinned crystals of YBCO show that the spin modu- 
lation occurs primarily along the a-direction, inferring that 
stripes are oriented along b [22]. This observation is en- 
tirely consistent with metallic screening channels being 
oriented along b. What remains unclear is the spacing of the 
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charge channels, which is more likely to be 2a from analy- 
sis of the phonon data [1], but assumed to be 4a by extrapo- 
lation of the diffraction measurements of pinned stripe 
materials to dynamic stripes [8,20]. 

The temperature dependence shown in Fig. 4 suggests 
that the gapped behavior of the Cu-0 bond- stretching mode 
is directly related to superconductivity. The changes in the 
intensity with temperature shown in Fig. 3 indicate that 
below T c even the A branch along the b direction develops a 
weak gap around 62 meV, so that the dispersion gap be- 
comes two-dimensional. This may be related to the two- 
dimensional nature of superconducting correlation [21]. It is 
important to note that the present result is consistent with 
the recent ARPES results on an untwined YBCO crystal 
[24] which show a strong anisotropy in the electronic struc- 
ture and a higher density of superconductive condensate in 
the direction perpendicular to the chain. These results sug- 
gest the possibility of superconductivity via the overscreen- 
ing mechanism as suggested by Tachiki and Takahashi 
[18], the phonon umklapp mechanism by Castro Neto [23], 
or the coupled spin-ladder mechanism of Sachdev [16]. 

For Lai. 85 Sro.i5Cu0 4 (LSCO) both the dispersion with a 
jump [1] as well as a continuous dispersion [25] were ob- 
served, suggesting that both of them exist in LSCO just as 
in YBCO [26]. Then the anisotropy in the screening may be 
a common feature of the cuprate superconductors, not a 
byproduct of the chain. In YBCO the orthorhombicity and 
the presence of the chain must favor the alignment of 
stripes in one direction throughout a twin domain, while in 
LSCO two types of domains with stripes along x and y must 
coexist. For large YBCO samples where twinning tends to 
obscure this anisotropy, the energy splitting of the bond- 
stretching phonons is fortuitous, and allows one to effec- 
tively observe a single domain lattice response for these 
strongly coupled modes. From the observations of the ani- 
sotropy of the lattice dynamics in this Letter and spin dy- 
namics and charge dynamics in the literature, it appears that 
all aspects of the low frequency dynamics are strongly ani- 
sotropic, as suggested by the stripe model. This places a 
stringent constraint to theories of high-temperature super- 
conductivity. 
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